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The Tolman-Oppenheimer-Volkov [TOV] equation constrains the internal structure of general 
relativistic static perfect fluid spheres. We develop several "solution generating" theorems for the 
TOV equation, whereby any given solution can be "deformed" into a new solution. Because the 
theorems we develop work directly in terms of the physical observables — pressure profile and density 
profile — it is relatively easy to check the density and pressure profiles for physical reasonableness. 
This work complements our previous article [Phys. Rev. D71 (2005) f 24307; gr-qc/0503007 wherein 
a similar "algorithmic" analysis of the general relativistic static perfect fluid sphere was presented in 
terms of the spacetime geometry — in the present analysis the pressure and density are primary and 
the spacetime geometry is secondary. In particular, our "deformed" solutions to the TOV equation 
are conveniently parameterized in terms of 5p c and 5p c , the finite shift in the central density and 
central pressure. We conclude by presenting a new physical and mathematical interpretation for the 
TOV equation — as an integrability condition on the density and pressure profiles. 
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I. INTRODUCTION 



The general relativistic static perfect fluid sphere has 
a long and venerable history that nevertheless continues 
to provide significant surprises [l|, 0, d, 0, d, H, 0, [1, H, 
[ToL I 111. [l2| . As emphasised in the review article by Del- 
gaty and Lake [l[ , while it is often relatively easy to write 
down explicit spacetime metrics that solve the differen- 
tial equation corresponding to a relativistic static per- 
fect fluid sphere, it is often much more difficult to check 
whether the corresponding pressure profile and density 
profile is "physically reasonable" — indeed there is often 
some confusion as to what the phrase "physically reason- 
able" might entail. 



This problem has if anything become even more acute 
with the recent introduction of "algorithmic" techniques 
that are capable of generating, and in principle classi- 
fying, all perfect fluid spheres 0, [H, H, uj- In view of 
this difficulty we have reformulated those "algorithmic" 
techniques in terms of several new "solution generating 
theorems" that apply directly to the TOV equation itself 
- and so directly give information about the pressure 
profile and density profile. Specifically, let us consider 
the well-known TOV system of equations, whose deriva- 



tion is now a common textbook exercise [III, : 

\p(r) + p(r)] [m(r) + 47rp(r) r 3 ] 
r 2 [l — 2m(r)/r] 



dp(r) 

dr 
dm(r) 
dr 



4-7T p(r) ', 



(1) 



(2) 
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We shall adopt units where G = c = 1, so that both m/r 
and p r 2 are "dimensionless" . 

Our basic strategy will be to assume that somehow we 
have obtained, or been given, some specific "physically 
reasonable" solution of the TOV equation in terms of 
the "initial" profiles po(r) and po(r), and then to ask how 
these initial profiles can be deformed while still contin- 
uing to satisfy the TOV equation. The output from the 
analysis will be several single-parameter, or sometimes 
multi-parameter, distorted profiles ps(r) and ps(r) that 
continue to solve the TOV equation for arbitrary values 
of the parameters E, along with some results regarding 
the "physical reasonableness" thereof. 

We shall use notation such as po,c and p c , and po,c 
and pc, to denote the central pressure, and central den- 
sity, before and after applying the distortion. Viewing 
the general relativistic static perfect fluid sphere as a 
first-approximation to a relativistic star, (though the dis- 
cussion could just as easily be applied to a perfect fluid 
planet) , we will concentrate on the region deep inside the 
stellar core, and specifically on the regularity conditions 
at the centre of the star, asking that the pressure and 
density remain positive and non-singular there. Observe 
that in seeking solutions of the TOV equation we are ef- 
fectively only demanding pressure isotropy — we make 
no claims as to constraints on, or even the existence of, 
an equation of state. 

We finally conclude by presenting a new physical and 
mathematical interpretation for the TOV equation — as 
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an integrability condition on the density and pressure 
profiles that allows us to explicitly write down the space- 
time metric for a relativistic perfect fluid sphere in terms 
of an iterated integral involving the density and the pres- 
sure profiles. 

II. LOCATING THE SURFACE 

As usual, we shall take the location of the first pressure 
zero p(r s ) = as defining the surface of the star, r s , 
though we should note a number of caveats: 

f. The "star" might not have a well-defined surface 
if p(r) ^ for all r — this corresponds to a "cos- 
mological" configuration where a localized star-like 
object is surrounded by an infinite diffuse "at- 
mosphere" which merges continuously into some 
asymptotic limit. While physically there is nothing 
necessarily wrong with such a configuration, they 
are certainly non-traditional [l[ . 

2. If the "star" is embedded in an asymptotically de 
Sitter or anti-de Sitter spacetime then the surface 
of the star is located by asking that the material 
pressure be zero. The total pressure at the sur- 
face is then determined by the cosmological con- 
stant, so the total pressure is nonzero at the sur- 
face: ptotai^s) = pa = —Pa ^ 0. Though at first 
this may seem a little alarming, at worst this sit- 
uation introduces some minor modifications of the 
standard analysis. 

3. The situation we shall implicitly be most interested 
in is where the spacetime is asymptotically flat and 
the surface occurs at some finite radius r s . To avoid 
needing too many special case technical analyses 
we shall always take our initial density and pres- 
sure profiles po( r ) and Po{ r ) to mathematically ex- 
tend to arbitrary r, including beyond r s . When we 
specifically want the physical density and pressure 
profiles we shall explicitly truncate them at the sur- 
face — for instance (in the pure vacuum case) by 



setting 

TO t runcated(r > r s) = m{r s ) = M, (3) 

and 

Ptruncatcd(r > r s ) = 0] (4) 



Ptruncatcd^ > T s ) =0. (5) 

These truncated profiles are certainly a solution 
to the TOV equation, and amount (in accordance 
with Birkhoff's theorem) to continuously joining 
the surface of the "star" to a portion of the exterior 
Schwarzschild geometry. 
4. If one wishes to embed a finite size star in de Sitter 
or anti-de Sitter spacetime then the relevant trun- 
cation conditions are 

47T 

m-truncated(r > r„) = m(r s ) + — p A (r - r s ), (6) 
and 

Ptruncatcd(r > T s ) = PA] (7) 



Ptruncatod(r > T s ) =PA = ~PA- (8) 

These truncated profiles are again certainly a solu- 
tion to the TOV equation, and amount (in accor- 
dance with a suitable generalization of Birkhoff's 
theorem) to continuously joining the surface of the 
"star" to a portion of the Schwarzschild-de Sitter 
[Kottler] geometry. 

We shall now use this general framework to develop 
several mathematical theorems regarding the TOV equa- 
tion, and subsequently study their physical implications. 



III. THE TOV EQUATION AS A RICCATI EQUATION 



Viewed as a differential equation for p(r), with m(r) [and hence p(r)] specified, the TOV equation is a specific 
example of a Riccati equation, for which the number of useful solution techniques is rather limited. (See for instance, 
Bender & Orszag [l5[ , or Polyanin and Zaitsev [l6| . The most important key results are collected in the appendix to 
this article.) Using these standard references it is relatively easy to show the following. 

Theorem (PI). Let po(r) and po(r) solve the TOV equation, and hold mo(r) — An J po{r) r 2 dr fixed. 
Define an auxiliary function go(r) by 

m (r) + A-K Po (r) r 3 
90{r)= r*[l-2mo(r)/r] " (9) 
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Then the general solution to the TOV equation is p(r) = po(r) + Sp(r) where 

Sp c 2m /r exp j-2 J .g drj 



8p(r) 



1 + Air8p c 



1 



exp < —2 / go dr > r dr 



(10) 



lo yjl - 2m Q /r 
and where Sp c is the shift in the central pressure. 

Proof. By using equation (|A2|) from the appendix we know that for any arbitrary real parameter k the quantity 



5p(r) — 



kexp 


rJ 


f m + Anpor 3 + &Trp r 3 ^' 




r 2 (l — 2mo/ r ) 


1 + ink j 


f r 


exp 


rJ 


fm + 4np r 3 + Sirpor 3 


dr 


1 — 2mo/r 


At 

r 2 (l - 2m /r) 



(11) 



is such that p(r) = po(r) + Sp(r) also solves the TOV equation. But by using an integration by parts one can simplify 
this to establish 



exp 



too + 47r po r 3 + 8irpo r 3 
r 2 (l — 2TOo/r) 



dr 



= yl — 2toq/V exp 



= y/l — 2mo/r exp 



m + 4irp r 3 
r 2 (l — 2TOo/r) 

-2 / ffo dr 



dr 



(12) 



Thus 



kyl — 2TOo/r 



exp 



6p(r) = 



2 / go dr 



l + ink 



2m /r 



exp < —2 / go dr > r dr 



(13) 



Finally, one fixes the constant k by first choosing the limits of integration to be to r, and then applying the boundary 
condition p(0) = p c at r = to deduce k = Sp c . Thus 



8p(r) — 



Sp c y/l- 2m /r exp | -2 J .g dr| 



1 + 4nSp c 



o y/i - 2m /r 



exp < — 2 / go dr > r dr 



(14) 



as claimed. 



□ 



A. Physically reasonable centre 

If the original po(r) and too (r) are "physically reason- 
able" , in particular if the central pressure po,c and den- 
sity p ,c are finite and positive, then m (r) = 0(r 3 ), so 
mo/r = 0(r 2 ), and <7o(r) = O(r). Thus all the integrals 
used above are finite and well-behaved (convergent) at 
the centre of the star. Consequently, the new central 
pressure will likewise be finite and well behaved; and will 
even be positive so long as Sp c > —po,c- So the central 
region of the new stellar configuration will automatically 
be "physically reasonable" . Since within the context of 
this theorem the density profile is unaffected, we also 
know that no horizon will be generated if we start from 
a horizon-free initial state. (Horizon formation is im- 



possible in the current situation: If it were to happen, 
it would be signaled by the existence of an r# such that 
2m(rij)/rH = 1, but by hypothesis p{r), and hence m(r), 
are held fixed. So if there is no horizon in the original 
configuration, one cannot generate a horizon by applying 
the current theorem.) 

Comment: With hindsight one can view this theorem PI 
as a consequence of theorem T2 of reference [4(. Some 
tedious algebraic manipulations will actually show that 
they are identical, up to a specific regularity choice. (In 
the notation of one must choose a ^ 0.) This regular- 
ity condition is now forced upon us because we want the 
pressure and density at the center of the star to remain 
finite. 

Equivalently, if the original po(r) and toq (r ) are "phys- 
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ically reasonable" , then (provided the relevant integrals 
are chosen to run from to r), the a = sub-case of 
theorem T2 of reference [3] yields a spacetime geometry 
for which g u — > at the origin, which we deem to not be 
physically reasonable. The key difference in the physics 
of the current article is that we now have a physics reason 
for restricting the solutions generated by theorem T2 of 
reference [ij , and now can thereby derive an explicit the- 
orem directly in terms of the shift in the pressure profile. 



B. Central redshift 

If one writes the spacetime geometry in the form 

dr 2 



ds 4 = -So(rydtf 



1 — 2m (r)/r 



+ r z d z n, (15) 



then the Einstein equations imply (see, for instance, 0, 

MM) 



Co(r) = exp <^ / g (r)dr 



(16) 



It is important to realise that with current conventions 
go(r) is defined by equation @, and that equation (JTC 
is a result of the Einstein equations. 
In particular, throughout the star 



CoW = Co, c exp | J g (r)dr }, (r < r s ), (17) 
and at the surface of the star 



Co{r s ) = Co,c exp \ J g {r)dr 



(18) 



where (o,c denotes the value of Co( r ) at the centre of the 
fluid sphere — which is effectively the central redshift of 
the fluid sphere. In fact for an idealized freely propagat- 
ing particle emitted from the centre 



Co,, 



1 



l + Z r 



(19) 



If one truncates the matter profile at r s , and and normal- 
izes using in the usual way using (truncated (°°) = lj then 
because of the Birkhoff theorem the geometry exterior to 
the star is Schwarzschild, and so one also has: 



Co(r s ) = y/l - 2M/r s 



(20) 



Consequently 



Co,c = Vl-2M/r s exp|- J g (r)drj . (21) 

Note that choosing the normalization (truncated (oo) = 1 
means that we are choosing the t coordinate to be proper 
time as measured at spatial infinity — which in general 
is well outside the surface r s of the star. For r < r s one 



could if one chooses rewrite equation (|17p in the com- 
pletely equivalent form 



(o(r) = y/l - 2M/r s exp |- J g Q (r)drj , (r < r s ). 

These considerations now allow one to rewrite the re- 
sult of theorem PI as: 



Sp(r) — 



5p c CL \/l-2m (r)/r 



Co(r) 



1 + An5p c C, 2 



0,c 



• dr 



/o v/l-2m /rCoW 2 

This has allowed us to eliminate the inner integral (in- 
volving g ) in terms of the g tt metric component. 



C. Shift in the surface 

From either equation (flU)) or ([2^]) . since all the other 
terms are positive, we see that the sign of Sp(r) is the 
same as the sign of 5p c , and so the same as the sign of 
the shift in the surface radius Sr s . This is most obvious 
for Sp c > where the positivity of Sp(r) implies that the 
surface definitely moves outwards. For small Sp c < the 
negativity of Sp(r) implies that the surface moves inwards 
- the tricky point is that for negative Sp c one could 
in principle (mathematically) encounter a pole at some 
finite r po ie, but this would then imply p(r po \ c ) = — oo, 
so the putative pole cannot occur inside the star — it is 
at worst a mathematical artifact located at r po i c > r s , 
outside the physical surface of the star. All in all, 



sign[(5r s ] = siga[Sp c ]. 



(24) 



We can say a little more if we restrict attention to 
infinitesimal shifts in the central pressure Sp c and take 
derivatives. The surface of the star is defined by 



p{p c ,r s (p c )) = 0, 



where 



Thus 



That is 



But now 
dp 

dp c 



p{Pc,r) 


= Po 


(r) + Sp(Sp c ,r). 


dp 




dp 


dr s 




dp c 


H 


dr 


r s ,p , c d P c 


PO,c 



= 0. 



dr s 
dp c 



dp/dp c \. 



dp/dr\ r 



d(Sp) 



r a ,Vo,c 



d(S Pc ) 



(25) 
(26) 

(27) 
(28) 
(29) 



r s ,8p c —0 



= y/l - 2M/r c exp I -2 J g drj ,(30) 
= Co,cA/l-2M/r 5 . (31) 
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Furthermore, from the TOV equation 

Po(r s )M 



dp 
dr 

r s ,po,e 
Combining, we see 



rf[l - 2M/r s 



(32) 



dr s 
dpc 



Co 2 c r 2 y/1 ~ 2M/r s 



above) the sign of 5r s is controlled, there is in general no 
a priori constraint on the sign of 3po(r s ) — po(r s ), so we 
cannot (without further assumptions) constrain the sign 
of the shift in compactness. This means that theorem PI 
by itself will not tell you if you are approaching or reced- 
ing from the Buchdahl-Bondi bound Xmaximum = 8/9. 
The best we can say is 



Po(r s ) M 

All terms appearing above are manifestly positive. 

D. Shift in the compactness 

If we now consider the "compactness" , defined as 
_ 2M _ 2m(r s ) 



(33) 



sign^x] = sign[3j0 O (r s ) - po(r 



sign[<Sr s 



(36) 



Note that the function {3po(r) — po{r)} — {3p(r) — p(r)} 
remains invariant under the conditions of theorem PI. 
Whether or not this quantity is positive or negative de- 
pends on one's a priori arbitrary choice of the initial so- 
lution profile {po(r), Po( r )} used as input into the theo- 
rem. If we now consider infinitesimal shifts in the central 
pressure we see 



(34) 



then, since m{r) = mo(r) is invariant under the current 
hypotheses, 

6 X=Y {3Po(r s ) - Po(r s )} r s 5r s + 0([8r s } 2 ), (35) 



dr s 3 
and via the chain rule 



{3po(r s ) - Po(r s )} r s , 



(37) 



where po(r s ) is the "average density" defined by 

3 



"^(^/("y^s). Unfortunately, while (as we have seen 



dx = 8tt {3p (r s ) - p (r s )} r s 3 C 2 e y/1 - 2M/r s 
dp c 3 



Po(r s ) M 



(38) 



J 



E. Shift in the spacetime geometry 



One might also reasonably ask what happens to the spacetime metric itself? As regards the radial g rr part of the 
spacetime metric, mo(r) is by construction fixed, and so g rr is unaffected. To calculate the "deformed" value of the 
the gtt — — C( r ) 2 component we note that from equation §§§ we have 



Sg(r) 



Ait Sp(r) r 
1 — 2mo(r)/r ' 



(39) 



whence 



5g(r) = 



47T 6p c r 
y/1 - 2m Q /r 



exp 



9o dr 



1 + Air8p c 



• dr 



o yl - 2m /r 



exp 



5o dr 



(40) 



That is 



r dr 



exp {-2 1 g dr 



2too/? 

Note that the Einstein equations imply, vide (fT6|) -([T7 |) . 

ln(C/Co) = ln(Cc/Co,c) + [ r Sgdr. 

Jo 

Integrating, 

f r rdr 

ln(C/Co) = ln(Cc/Co,c) + In <^ 1 + Ak8 Pc \ ^ = , exp 



o \J\ - 2m /r 



9o dr 



(41) 



(42) 



(43) 
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The only tricky part of the computation lies in getting the overall normalization correct. Since go (r) is defined by 
([5]) directly in terms of mass and pressure profiles, the above formula provides a completely explicit formula for the 
(multiplicative) shift in £, and so for the (multiplicative) shift in g u - Explicitly, 

«"-«"»L wt ,a,r „ : d \ 1. (44, 



Cc Co,c [ ' c Jo \A~ 2m /r £ (r) 

In view of these formulae for £(r) we can rewrite the shift in pressure in the more compact manner 



s ( \ Sp c y/l- 2m Q r Co, c Cc , AK s 
Co(r) CM 

Though significantly more compact, the trade-off in this formula is that £(r) and £o( r ) are themselves quite complicated 
functions of the mass and pressure profiles. 

Normalization of theorem PI versus theorem T2: With the above normalization results now in hand, we see that 
theorem PI corresponds to theorem T2 of [4j with the integrals running from to r, with a — > Cc/Co,o while the 
other parameter e occurring in T2 can be physically identified by taking e/er — > 47r 5p c Co c - 



F. Two additional theorems 

In the specific case considered in this section, where m(r) is held fixed, one can obtain two additional related 
theorems by using equations (|A3j) and ()A4j) of the appendix. Specifically: 

Theorem (Plb). For a fixed specification of po{r), letp\{r) andp2{r) be two distinct solutions of the TOV equation. 
Then the general solution (for X an arbitrary real parameter) is 



X exp 



-^p—r- rdr\pi(r) + (1 - A) cxp { - f 4 ^ P2 r dr\ p 2 (r) 



p(r) = ^ JIJ 1 r 2m0/ : I >- Jl ' r 1 — " u/ : ^ . (46) 

Aexp {"i T^A :rdr } + (1 - A)exp {-i T^oTr rdr } 

The central pressure is 

Pc = A p c ,i + (1 - A) p C)2 . (47) 

Note that this theorem no longer requires any nested integrations, at the cost of needing two specific solutions as 
input. A minor technical issue is that pi(r) and p-z(f) have to solve exactly the same Riccati equation over the entire 
domain of interest — this leads to technical problems if we truncate the density and pressure profiles at the surface of 
the star. We will either be forced to work on the intersection of the truncation regions, that is r S [0, min{r S! i, r S) 2}], 
or work with un-truncated pressure profiles pi(r) and P2(r) that have been extended into the region where pressure 
is permitted to become negative. 

Theorem (Pic). For a fixed specification of po(r), letp\(r), pi(r), and ps(r) be three distinct solutions of the TOV 
equation. Then the general solution (for X an arbitrary real parameter) is 

( \ = Api(r) [p 3 (r) -p 2 {r)\ + (1 - A) p 2 (r) [p 3 (r) - Pi{r)\ 
PK) A[p 3 (r)-p 2 (r)] + (l-A) \p 3 (r) - Pl (r)] ' 1 ' 

The central pressure is then 

p = Ap c ,l [Pc,3-Pc,2] + (1 - A) p e , 2 [Pc,3 -Pc,i] / 4g s 
A [Pc,3 ~Pc,2] + (1 - A) \p c<3 -Pc,l] 

Note that this last theorem no longer requires any integrations whatsoever, though now at the cost of needing three 
specific solutions as input. As in the preceding theorem, p\(r), pi(r), and pa(r) have to solve exactly the same 
Riccati equation over the entire domain of interest — for stars with a well defined surface we will either be forced 
to work on the intersection of the truncation regions, that is r e [0, min{r Si i, r s ^ 2 , ^,3}], or work with un-truncated 
pressure profiles Pi(r), p 2 {r), and ps(r) that have been extended into the region where pressure is permitted to become 
negative. 



IV. CORRELATED CHANGES IN DENSITY sure profiles. 

AND PRESSURE 



The second main theorem we develop can be obtained 
by looking for correlated changes in the mass and pres- 
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Theorem (P2). Let po(r) and po(r) solve the TOV 
equation, and hold go fixed, in the sense that 

_ m (r) + Aitpo{r)r 3 _ m(r) + 4irp(r)r 3 
9o[r > ~ r 2 [l-2m (r)/r] ~ r 2 [l - 2m(r)/r] ' [ ' 

Then the general solution to the TOV equation is given 
by p(r) =po(r) + 5p(r) and m{r) = mo(r)+Sm(r) where 

Airr 3 5p c f f 1 - r g \ 
° m ( r ) = 7rn~, 19 ex P i 2 / 9o z — ; dr > , (51) 



3 [l + rgo] 2 



8p(r) 



o 1 + rgo 



5m 1 + 8irpor 2 



47rr 3 1 — 2mo/r 



(52) 



Here 5p c is the shift in the central density. 
Explicitly combining these formulae we have 



Sp(r) 



6p c 1 + 8irp r 2 



[1 + r go} 2 1 - 2mo/r 

x exp { 2 / .go ] ; * 9 ° dr 
{ Jo l + rg 



Proof. Since 



m(r) — 



go(r)r 2 — 4irp(r)r 3 
1 + 2rg a (r) 



(53) 



(54) 



we have in particular 



An6p(r) r 3 3 
6m(r) = — L —- = -<iiT8p{r)r A 



1 + 2rg (r) 



1 — 2mo/r 
1 + 8irpor 2 



so that 



Sp(r) 



8p(r)r 3 



1 d 

r 2 dr V 1 + 2r , go('') 



(55) 



(56) 



Now consider the TOV equation, or more precisely, the 
change in the TOV equation: 



dSp{r) 
dr 



= ~[Sp(r) +8p(r)} g (r). 



(57) 



Combining these last two differential equations yields a 
linear homogeneous differential equation for 8p(r), which 
is easily integrated. The quoted form of the theorem 
results after an integration by parts. □ 



A. Physically reasonable centre 

If the original density and pressure profiles are "physi- 
cally reasonable" , then likewise the new density and pres- 
sure profiles will be well behaved — at least for a finite 
region including the origin. In particular, since regular- 
ity of the original solution implies go{r) — 0(r) at the 
origin, the only integral we have to do for this theorem 
is guaranteed to converge at the lower limit r = 0. So 



the central region of the new stellar configuration will 
automatically be "physically reasonable" . Note that for 
theorem P2 



s Pc = - 5 -^. 



(58) 



If the central pressure increases 5p c > 0, then for the sit- 
uation envisaged in this theorem the mass profile always 
decreases, implying that an event horizon never forms. 

Comment: With hindsight one can also view this theo- 
rem P2 as a consequence of theorem Tl of reference 0]. 
The key difference now is that we have an explicit state- 
ment directly in terms of the shift in the pressure profile. 
To see this, write go — Co/Co , substitute into Tl and 
apply boundary conditions to deduce A — > 87r 5p c /3. 



B. Shift in the surface 

Note that under the conditions of theorem P2 we have 

sign[(5p(r)] = — sign[<5m(r)] = — sign[#p c ]. (59) 

So if the central density goes up, the pressure everywhere 
goes down, and the radius of the surface of the star (being 
defined by the first pressure zero) must decrease. That 
is 



sign[<5r 



-sign[<5p c ] = +sign[<5p c 



(60) 



If we now consider infinitesimal changes in central pres- 
sure we see that equation (|28|) is unaltered. Furthermore, 
under the hypotheses of theorem P2 we have 



dp 
dp c 



d(6 P ) 



r a ,Po,a 



d(8 Pc ) 



r s ,Sp c — 

1 



[1 + r s go{r s )] 2 [1 - 2M/r 
1 - rg 



x exp 



9o 



1 



dr 



rg 



(61) 



(62) 



which is manifestly positive. Thanks to the TOV equa- 
tion we again have equation (|32[) , so combining the above 
leads to 

dr s 
dp c 



[1 + r s go(r s )} 2 po(r s ) M 
x exp i 2 



1-rgo 
9o dr 



1 



rgo 



(63) 



All terms appearing above are manifestly positive, as was 
the case for theorem PI. 



C. Shift in the compactness 

The compactness x — 2?n(r s )/r s now changes both due 
to changes in m(r) and due to changes in r s . For small 



8 



but finite changes we have 
2 Sm(r s ) 8tt 



SX = 



+ IT {3po(r s ) - Po(r s )} r s Sr s + 
r 6 



It is more efficient to re-write the compactness as 

2m(p c ,r s (p c )) 

XKPc) = f-^ 

rs(Pc) 

and then differentiate. We see 
2 dm(p c ,r) 



dp c 



dp c 



(64) 
(65) 

(66) 



. 8vr _ dr s {p c ) 

+ Y i 3 Povs) - Po{r s )} r s — ^— — 



PO.c 



But under the hypotheses of theorem P2 we have 



dm(p c ,r s 



dp c 



[l + r s g f 



(67) 



x exp { 2 / g - T 90 dr 
o l + rgo 



Now, combining this with the results of the previous sub- 
section 



dx 
dp c 



M\l + r s g ] 



exp <2 g Q — dr 



o ^ + rg 



r s 3p {r s ) 



(68) 



While the first few terms are guaranteed positive, there 
is no such guarantee as to the sign of the last term. Thus 
without further hypotheses we again cannot constrain the 
sign of the shift in the compactness \ — 2m,(r s )jr a . 



V. COMBINING THE PREVIOUS THEOREMS 



One can obtain additional and more complicated the- 
orems by iteratively applying theorems PI and P2, in 
a manner similar to the discussion of reference [H . This 
iteration will yield 2-parameter generalizations of theo- 
rems PI and P2. 



Theorem (P3). Let po(r) and po(r) solve the TOV 
equation. Apply theorem PI followed by theorem P2. 
Let us define intermediate quantities 



5pi(r) 



Ap y/l- 2m /r exp j— 2 J g dr j 



1 + 47T Ap 



rdr f r ^ 

o y/l - 2m /r 1 Jo 

(72) 



9i{r) = go(r) + 



Ait 5pi(r) 
1 — 2mo/r 



(73) 



Then p(r) = po(r) +6p(r) and m(r) — mo(r) + 5m(r) are 
also solutions of the TOV equation, where 



Mr) = -p (2 IZIML dr \ , (74) 



3[l + r S i] s 



and 



5p(r) = Spi(r) 



1 + rgi 



8m 1 + 87r[po + 5p\]i 



Aur 6 



1 — 2mo/r 



(75) 



D. Shift in the spacetime geometry 

What happens to the spacetime geometry under the 
application of theorem P2? By construction the quantity 
g(r) = go(r) is unaffected, and therefore 



CO") Co(r) 



Cc 



Co.. 



(69) 



It is only the g rr component of the spacetime metric that 
changes significantly: 



1 



g r 



[gojrr 1 - S[2m(r)/r] ' 



(70) 



where we already know 



2m(r) 



8nr 2 Sp c f f r l-rg 

—-— exp ^ 2 / g — dr 

3 [l + r. 9o ] 2 I Jo 1 + ^50 



Here 8p c is the shift in the central density, and the total 
shift in the central pressure is given by 



5p c = Ap - 



(76) 



The two parameters Sp c and Ap can be specified indepen- 
dently. 

Theorem (P4). Let po(r) and po(r) solve the TOV 
equation. Apply theorem P2 followed by theorem PI. 
Then p(r) = po(r) + Sp(r) and m(r) = ma(r) + 6m(r) 
are also solutions of the TOV equation, where 



47rr 3 Sp c , 
om[r) = — exp < 2 



3 [1 + r 50] 



fgo'—^drl (77) 
'0 l + ^5o 



(71) and 



Sm 1 + Sirpor 2 



4nr 3 1 — 2rrio/r 



Ap Co 2 c V 1 "2[ m o + '5m]/r 



1 + 4ttAp Co,, 



r 



r dr 



o ^/l-2[m + <5m]ACo(r) 



(78) 



-ffere 8p c is the shift in the central density, and the total 
shift in the central pressure is given by 



Sp c = Ap - 



(79) 



The two parameters Sp c and Ap can be specified indepen- 
dently. 

Comment: Note that these two theorems provide two- 
parameter generalizations of the TOV solution (p , Pa) 
that one starts from. These theorems are closely related 
to theorems T4 and T3 of reference [H . Technically they 
are equivalent to the case of theorems T4 and T3 

of reference [3] , where a ^ is now a physical restriction 
we place on the solution by demanding that the density 
and pressure be well behaved at the centre of the fluid 
sphere. 



VI. THE TOV EQUATION AS AN ABEL 
EQUATION 

If, on the other hand, we think of the pressure profile 
p(r) as fixed, than we can rearrange the TOV system of 
equations as a differential equation for m(r), specifically: 



dm(r) 
dr* 



-^(r^- ^l-^/r] ^ (80) 
m[r) + 47rp(r)r' r> dr 



This is now an Abel equation (2nd type, class A). 

If one were able to develop a solution generating the- 
orem based on this equation it would in many ways be 
the most natural companion to theorem PI. Unfortu- 
nately, despite numerous attempts at developing solu- 
tion generating theorems based on this observation, we 
have no progress to report. This is ultimately due to the 
fact that mathematically Abel equations are considerably 
more difficult to deal with than Riccati equations. 



VII. THE TOV EQUATION AS AN 
INTEGRABILITY CONDITION 

We now develop an alternative point of view regarding 
the physical and mathematical significance of the TOV 
equation: In any spherically symmetric spacetime there 
are formally three algebraically independent components 
of the Einstein equations. Let us now formally solve two 
of the Einstein equations by building them into the cho- 
sen form of the line element, and use the remaining Ein- 
stein equation to derive (in the case of perfect fluids) an 



integrability condition — this provides us with a new and 
different way of interpreting the TOV equation. 

Specifically, let us now consider two arbitrarily speci- 
fied functions p(r) and p r (r), and define 



m(r) 



Anp{r)r 2 dr. 



(81) 



Now consider the line element 



ds 2 



exp 



4-7r[p(r) + p r (r)]r 
1 — 2m(r)/r 



dr 



1 - 
dr 2 



2m(r) 



dt 2 



1 — 2m(r)/r 



r 2 (d0 2 



>d</> 2 ). (82) 



That is, the line element is specified completely by a 
double iterated integral involving the arbitrary functions 
{p(r),p r (r)}. Then two of the Einstein equations are 
trivial, they are "built in" to the chosen line element: 



Gq = 8ttGn p\ 
Gff = 8itGn p r . 



(83) 
(84) 



These two equations permit us to physically interpret 
p(r) as the density and p r (r) as the radial pressure. 

The sole remaining Einstein equation depends on the 
Ggg = G ^ component of the Einstein tensor. Now, it is 
a purely geometrical statement that 



Goo — Gj 



(85) 



Gff + 47rr 



dp r {p + Pr)(m + 4irp r r 3 ) 
dr r 2 (l — 2m/r) 



If we enforce isotropy, Ggg = G^ — Gff, then without 
further ado we deduce the TOV equation 



dp r (p + p r )(m + 4np r r 3 ) 
dr r 2 (l — 2m/r) 



(86) 



which we now see can be interpreted as the integrability 
condition allowing us to solve the two equations (|83[) - 
([84)1 for a perfect fluid source. Of course once we have 
the TOV equation, however derived, we can apply the 
Riccati equation analysis developed earlier in this article. 

A bonus of setting things up this way is that it is im- 
mediately clear what happens for anisotropic fluids. The 
third Einstein equation now yields 



Pt = Pr 



dpr (p + Pr)(m + 4:Trp r r 3 
dr r 2 (l — 2m/ r) 



(87) 
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This can either be viewed as a way of calculating the 
transverse pressure p t in an anisotropic fluid, or can be re- 
arranged to yield the well-known anisotropic TOV equa- 
tion 

dpr = (p + p)(m + 4irpr 3 ) 2(p t - p r ) 
dr r 2 (l — 2m/ V) r 

VIII. DISCUSSION 

The purpose of this article has been to develop several 
"physically clean" solution generating theorems for the 
TOV equation — where by "physically clean" we mean 
that it is relatively easy to understand what happens to 
the pressure and density profiles, especially in the vicinity 
of the stellar core. In many ways this article serves as a 
companion paper to reference [J], where related results 
were derived in terms of the spacetime geometry. These 
new results refine our previous results, in the sense that 
it is now clearer when a mathematical solution of the 
isotropy equations are "physically reasonable", at least 
in the region of the stellar core. 

Indeed one important message is that in theorems T2, 
T3, and T4 of reference [4[ one should take a ^ if one 
wishes the pressure and density to remain well behaved 
at the center the of the fluid sphere. 



Reference [4J also contains a number of other inter- 
esting results regarding the use of "solution generating" 
theorems to classify perfect fluid spheres — sometimes 
new solutions are obtained, sometimes old solutions are 
recovered in a new context, and theorems T1-T4 can 
be used to generate a web of interconnections between 
known and new perfect fluid spheres. These comments 
also apply mutatis mutandis to the present theorems, and 
we direct interested readers to reference [i[ for further de- 
tails. 

In closing, we reiterate that the general relativistic 
static perfect fluid sphere, despite its venerable history, 
continues to provide interesting surprises. 
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APPENDIX A: KEY RESULTS ON RICCATI EQUATIONS 

The general Riccati equation is 

a(r) + P(r) p(r) + j(r) p(r) , (Al) 



dr 

for specified functions a(r), (3(r), 7 (r). While no completely general solution from first principles exists, if we are 
given one specific solution po(j) then the one-parameter general solution (parameterized by an arbitrary real number 
k) may be written 



k exp j 




f[2 7 (r)po(r) + /3(r)]drj 


\ 


1 — k J j(r) exp | 


U 


/"[2 7 (r)p (r)+/3(r)]drj 


|>dr 



P{r) = Po (r) + ^-r- J —, • ( A2) 



This can easily be derived, for instance, from the argument sketched in the standard reference Bender & Orz ag [la |. 
alternatively an equivalent explicit statement be found in the reference handbook by Polyanin and Zaitsev [la ]. A 
second useful result is that if we know two specific solutions pi(r) and P2(r) then the one-parameter general solution 
(now parameterized by an arbitrary real number A) may be written [161 ]: 



A exp | 


^ -J 7(r)pi(r)dr 


>Pi(r) + (l-A) expj 


n 


f 7 (r)p 2 (^)drj 




A 


exp<j 


{'J 


f 7(r)pi(r)drj 


[ + (1 - A) exp < 




^ j{r)p2(r)dr 


\ 



P( r ) = — — 7 : 1 s r^-r ■ ( A3 ) 
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Finally, if we know three specific solutions pi(r), P2{ r ), and p% (r) then the one-parameter general solution (again 
parameterized by an arbitrary real number A) may be written [16| : 

= Api(r) [p 3 (r) -p 2 {r)\ + (1 - A)p 2 Q) [p 3 (Q ~ Pi(r)} 
PU A [p 3 (r)-p 2 W] + (l-A)tp s (r)-p 1 (r)] ' 1 j 

Note that these three forms of the general solution require two, one, or zero integrations respectively, at the cost of 
needing one, two, or three pre-specified input solutions. 
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